Quasi-universal finite-T scaling in gapped one-dimensional quantum magnets. 
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Temperature dependencies of gap energies and magnon lifetimes are measured in the quasi- 1- 
dimensional 5=1/2 gapped quantum magnets IPA-CuCla and SUI-CU2CI4 using inelastic neutron 
scattering. The results are compared to those found in literature for 5* = 1 Haldane spin chain 
materials and to theoretical calculations for the 0(3)- and 0{N)- quantum non-linear a-models. It 
is found that when the T = energy gap A is used as the temperature scale, all experimental and 
theoretical curves are identical to within system-dependent but temperature-independent scahng 
factors of the order of unity. This quasi-universality extends over a surprising broad T range, at 
least up to kT 1.5A. 
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One-dimensional (ID) physics is unique in that colli- 
sions between counter-propagating particles becomes un- 
avoidable, regardless of their mutual interaction strength. 
This circumstance has drastic consequences, such as the 
collapse of the Landau Fermi Liquid in ID conductors 
and of long range order in gapless spin chains [1)]. In 
these extreme examples the quasiparticle picture totally 
breaks down. Not so in gapped quantum spin liquids, such 
as Haldane spin chains [2] or spin ladders 0, [3]. Here 
the quasiparticles (magnons) persist in the limit T — > 0, 
since mutual collisions are rare due to their exponentially 
small density. Interactions do become important at ele- 
vated temperatures thanks to a thermal excitation of a 
large number of magnons. The result is a reduction of 
magnon lifetimes and a renormalization of their energies. 

In the simplest case the energy gap A is large compared 
to the dispersion bandwidth and magnons are localized. 
The finite-T effects are then well accounted for by the 
Random Phase Approximation (RPA) [^, @]. This ap- 
proach ignores correlations and breaks down if the quasi- 
particle are highly mobile, with a velocity v ^ A. For- 
tunately, the latter regime is open to investigation us- 
ing field-theoretical methods. Moreover, as mentioned 
above, in ID the details of the interaction potential be- 
come irrelevant. For kT ^ A a universal theory with 
a single energy scale A can be derived from the large-5 
mapping of the Heisenberg Hamiltonian onto the 0(3) 
non-linear sigma- model (NLSM) 0, 0]. Until recently, 
these predictions and the domain of their applicability re- 
mained largely untested experimentally. Of the strongly 
ID gapped quantum magnets, finite-T data existed only 
for 5=1 Haldane spin chains [9, l3, [HI, 0, E^, where 
A = v/A ~ 6 [14]. In the present work we explore the 
temperature dependence of inelastic neutron spectra in 
the 5=1/2 spin ladder system IPA-CuCla (A - 2) 
ll [H, and the 4-leg spin tube Sul-CuaCU (A > 25) 
18| . We find a striking quasi-universal renormalization 
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FIG. 1: (a) Typical background-subtracted constant-g scans 
collected at the ID zone-center (0.5, 0, 0) in the spin ladder 
compound IPA-CuCla {v/A — 1.9) at different temperatures 
(symbols). Solid lines are fits to the data, as described in 
the text, (b) False-color representation of inelastic intensity 
measured at the ID AF zone-center in IPA-CuCla plotted as 
a function of temperature and energy transfer. Symbols are 
the measured temperature dependence of the gap energy. The 
solid line is a guide for the eye. 



of the magnon gaps and lifetimes that persists to surpris- 
ingly high temperatures. 

As described in detail in [l^, the magnetic properties 
of IPA-CuClg are due to 5 = 1/2 ladders of Cn'^+ ions 
that run along the crystallographic a axis of the PI crys- 
tal structure. Exchange interactions along the ladder legs 
are antiferromagnetic (AF) . Pairs of spins on each ladder 
rung are correlated ferromagnetically. The ID AF zone 
center is at h — 0.5, where (h,k,l) are components of a 
wave vector q. Due to residual interactions between lad- 
ders along the c axis, at low temperatures the spin gap 
at ft, = 0.5 varies between A = 1.17 meV for ^ = to 
A = 1.8 meV at I = 0.5. There is no detectable disper- 
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FIG. 2: Same as Fig. 1, but for the (0,0,0.5) ID 4-leg spin- 
tube material SUI-CU2CI4 {v/A > 20). 



sion of magnetic excitations along the b axis. Treating 
the weak inter-ladder coupling at the RPA level, one can 
estimate the intrinsic ID spin ladder gap Aq = 1.5 meV. 
The magnon dispersion along the ladder axis is parabolic 
at the ID AF zone-center and characterized by a veloc- 
2.9 meV, as defined by Eq. 2 in 



ity V — 2.9 meV, as defined by Eq. 2 in [17|, so that 
X — v/Aq — 1.9. The other Cu^"*" based material that 
we discuss here is SUI-CU2CI4 The corresponding 

topology of its magnetic interactions is rather complex, 
but comes down to that of 4-leg spin tubes with very 
strong AF leg coupling along the c axis. Interactions be- 
tween adjacent spin tubes are negligible and the system 
is exceptionally ID. Due to geometric frustration, the ID 
AF zone-center is located at slightly incommensurate po- 
sition / = 0.5 — C, where C = 0.022(2). ( amounts to only 
about 20% of our experimental wave vector resolution, 
and can be safely ignored in the context of the present 
study. At low temperatures the gap in SUI-CU2CI4 is 
Aq = 0.55 meV. Dispersion along the c* direction is very 
steep, with v = meV. Thus, for SUI-CU2CI4, A = 25 
is 12 times larger than for IPA-CuCla and 4 times larger 
than for Haldane chains. 

Temperature-dependent inelastic measurements were 
performed using 3-axis neutron spectrometers. The data 
on IPA-CUCI3 at momentum transfer q = (0.5, 0, 0) were 
collected on the IN22 instrument at ILL using 5 meV 
fixed-final neutron energy, a pyrolitic graphite (PG) 
monochromator, a horizontally focused PG analyzer, and 
a Be filter after the sample. The background was mea- 
sured independently at each temperature away from the 
ID AF zone-center, at q = (0.3, 0, 0) and q = (0.7, 0, 0). 
The sample consisted of 20 co-aligned fully deuterated 
single crystals of total mass 3 g. Sample environment 
was a standard He-flow cryostat. A similar configura- 
tion on the NG-5 SPINS instrument at NIST was used 
with 3.7 meV fixed final energy neutrons to collect data 
around (0.5, 0, 0.5) in IPA-GuGlg and at (0, 0, 0.5) in Sul- 
CU2CI4. For the latter material, we utilized a ~ 2 g as- 



sembly of 15 co-aligned deuterated single crystals and a 
^He-'^He dilution refrigerator. 

Typical background-subtracted energy scans collected 
at the ID AF zone-centers in IPA-CuCls and Sul- 
GU2CI4 at different temperatures are shown in Fig. [1^ 
and Fig. [2^, respectively. Series of such scans are assem- 
bled into false-color plots in Figs[T]3 and Fig.Oa. Here the 
measured intensity at each temperature was normalized 
by its peak value. The observed gap modes broaden with 
temperature, and the gap energies increase. The data 
were fit to a model cross section function that was nu- 
merically convoluted with the resolution function of the 
spectrometer. The dynamic structure factor was written 
in the two-Lorentzian form, as given by Eqs. 6 and 7 
in 10]. Lorentzian line shapes are consistent with the 
theoretical results of Q . The excitation width P was as- 
sumed to be q-independent. The dispersion relation was 
written as [fixjo(q)]^ = A^ + i;^(qd — tt)^, with d = a for 
IPA-GuGla and d = c for SUI-CU2GI4, respectively. The 
parameters P, A, and an overall intensity prefactor were 
refined to best fit the scans measured at each tempera- 
ture. The velocities v were fixed at their previously deter- 
mined low-temperature values. Typical fits are plotted 
as heavy solid curves in Figs. [1^ and [2^. The tempera- 
ture dependencies of the gap energy A(r) obtained from 
the fits are plotted in symbols in Figs. [It and [2)3. In 
Fig.[3]we plot the dimensionless variation of A, defined in 
quadrature as S(T) = v^A2(T) - A2(0)/Ao, to empha- 
size the low-temperature region. The abscissa axis shows 
the reduced temperature r = kT/Aq, k being the Boltz- 
man's constant. The r-depcndcnce of the relative excita- 
tion half-width -f{T) = P(r)/Ao is shown in Fig. H [13]. 
Note that for IPA-CuCls we obtain a perfect overlap be- 
tween measurements at the minimum and maximum of 
the transverse dispersion, validating our expectation that 
the temperature dependence is an intrinsic ID effect. 

Central to this study is a comparison of our results for 
SUI-CU2GI4 and IPA-CuCls to those reported in fitera- 
ture for other ID quantum spin liquids. Extensive finite- 
T data exist for the S = 1/2 dimer systems Cu(N03)2 
21 1 and TlCuCla . Unfortunately, in the former com- 
pound V <^ A, while the latter is by no measure one- 
dimensional. Among the known materials that do fit our 
requirements, neutron data exist only for S = 1 Haldane 
spin chain systems. CsNiCla is one of the oldest known 



prototypes |23l . |24| . Inter-chain interactions in this com- 



pound are significant, and it actually orders in 3D at T = 
4.8 K. At the 3D zone-center the gap softens at the transi- 
tion point [10!], but one can hope to retrieve ID behavior 
at certain special wave vectors 12|. The intrinsic ID gap 
energy in GsNiGla is A(0) ~ 0.9 meV. In Figs, ^and 
|4] we plot the temperature-dependent data from [l2l| in 
star symbols. The organo-metallic complexes NENP [25| 
and NINAZ ll| are excellent ID systems, but they are 
affected by the rather strong magnetic anisotropy. This 
anisotropy splits the triplet of Haldane gap modes into 
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FIG. 3: Reduced gap energy 5 plotted against reduced 
temperature r for diverse quasi ID quantum antiferromag- 
netic materials (symbols) and analytical models (lines). In- 
set: same data scaled to show universal behavior. CsNiCls, 
NENRNINAZ and Y2BasNi05 data are from Refs. [H, H [n| 
and[l^, respectively. Calculations for the 0(3) and 0{N) non- 
linear (T-models are from Refs. [l^ and[20l. correspondingly. 



a low energy transverse-polarized doublet and a longitu- 
dinal higher-energy singlet 26:j|, with gaps A_l and A||, 
respectively. As argued in [19|], in this case one can still 
make a meaningful comparison with Heisenberg systems 
by focusing on the two lower modes. For NENP and NI- 
NAZ A_L = 1.1 meV and = 3.6 mcV, respectively. In 
Figs. El and H we plot S{T) ^ ^A2^(r) - A2_(0)/A_l(0) 
and 7(T) = T^{T)/A^{0) vs. r = r/Ax(0), with up 
and down triangles for NENP [§] and NINAZ fllj. re- 
spectively. The most recent and extensive study [13| is on 
Y2BaNi05 that is known to be an excellent ID Haldane 
gap system with A = 8.6 meV and only weak anisotropy 
[2 71 ] . In Figs. [3] and [4] the corresponding data from [13| 
are plotted in square symbols. As mentioned above, for 
S = 1 Haldane spin chains numerical calculations predict 
A = 6.2 [ij], but it is also useful to quote the experimen- 
tall ratio A = v/A = 6.4 for CsNiClg pi, A = 6.1 for 
NENP and A = 8.1 for YaBaNiOs H^. 

It is instructive to compare the experimental data 
on IPA-CuCls, SUI-CU2CI4 and the other materials de- 
scribed above to field-theoretical results. The solid line 
in Fig. [3] shows a one- loop calculation of the temperature 
dependence of the gap energy in the 0(3)-NLSM from 
[l9|. For this model thermal broadening of excitations 
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FIG. 4: Reduced energy width 7 of the gap excitations plot- 
ted against reduced temperature r for diverse quasi ID quan- 
tum antiferromagnetic materials (symbols, as in Fig. [3]) and 
the 0(3) non-linear u-model (dash-dot line, Ref. @). 



was calculated in [sj , and is plotted in a dash-dot line in 
Fig. m More results are available for the large- limit 
of the NLSM that is more easily tractable. 0(A^)-NLSM 
calculations for A(T) were reported in [2o!| and are plot- 
ted in a dashed line in Fig. [3] 

The thermal renormalization of gap energy remains 
within a factor of two for all the diverse materials and 
models considered. From the theoretical viewpoint, the 
increase is driven by a repulsion between thermally ex- 
cited quasiparticles 0, S [I^. It is proportional to 
the correlation length A and the quasiparticle density 
p. At low temperatures the latter universally scales as 
p oc -y/r/Aexp (— 1/r) for all ID systems with a dynam- 
ically generated gap. As a result, for kT A one can 
expect all measured S{t) dependencies, as well as the 
ones calculated for the 0(3)- and 0(A^)-NLSM, to be the 
same, to within a system-dependent but temperature- 
independent scaling factor a. The latter will reflect the 
details of quasiparticle interactions in each particular 
model. Experimentally, this is indeed the case in the 
entire studied temperature range. The inset in Fig. [3] 
shows such a scaling plot, where we have arbitrarily se- 
lected a = 1 for the 0(A^)-NLSM, a = 1.3 for the spin 
ladders in IPA-CuClg, a = 1.35 for Sul-CuaCU, and 
a = 2.35 for all the Haldane spin chains. The data on 
CsNiCla were excluded from the scaling, as they are in- 
consistent with those on other S = 1 chains, due to the 
strong 3D interactions in that material [lo| . 

A similar quasi-universal scaling can be also expected 
for the excitation widths. The latter is given by the in- 
verse average time between collisions among thermally 
excited quasiparticles and is thus proportional to p and 
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thermal average of quasiparticle velocity v oc v^T/A Q. 
The experimental curves in Fig. |4] are indeed very close 
for all systems under consideration, and a near-perfect 
collapse (inset) is obtained with scaling factors /? = 1 for 
the 0(Ar)-NLSM, f3 = 1.2 for the spin ladders in IPA- 
CuCla, (3 = 0.8 for Sul-CuzCU, and l3 = 1.6 for Haldane 
chains. 

Deviations from unity of the factors a and /? reflect the 
inadequacy of the NLSM mapping. Indeed, it is rigorous 
only in the limit S » 1. As evident from recent numeri- 
cal studies [28| , that model has only limited quantitative 
applicability to S — 1/2 and S — 1 chains. An alterna- 
tive description of the 5 = 1/2 ladder contains quasipar- 
ticle interactions as a free parameter [29| . and may be 
more versatile. Of course, any quasi-universal behavior 
with the energy scale set by A will have a solid theoret- 
ical justification only in the limit kT ^ A. A different 
universality is expected for v 3> kT 3> A, where the 
field theoretical description is still justified, but the spin 
gap can be considered negligible. Here the only remain- 
ing energy scale is the temperature itself and x"(7r,a;) 
is a function of cu/T. Most of our data actually lie in 
the crossover region nT ~ A, while access to the limit 
kT ^ a is hindered by the limited experimental resolu- 
tion. The two central experimental findings are thus the 
quasi-universal T/ A scaling in vastly diverse models and 
real small-S* materials, and its persistence well beyond 
the low-temperature limit, up to at least kT ^ 1.5 A. 
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